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The spin-1 anisotropic antiferromagnet NiCl2-4SC(NH2)2 exhibits a field-induced quantum phase 
transition that is formally analogous to Bose-Einstein condensation. Here we present results of 
systematic high-field electron spin resonance (ESR) experimental and theoretical studies of this 
compound with a special emphasis on single-ion two-magnon bound states. In order to clarify some 
remaining discrepancies between theory and experiment, the frequency-field dependence of magnetic 
excitations in this material is reanalyzed. In particular, a more comprehensive interpretation of the 
experimental signature of single-ion two-magnon bound states is shown to be fully consistent with 
theoretical results. We also clarify the structure of the ESR spectrum in the so-called intermediate 
phase. 



PACS numbers: 

I. INTRODUCTION 

The organic compound NiCl2-4SC(NH2)2 (known as 
DTN) is a gapped quasi-one-dimensional antiferromag- 
net with easy-plane anisotropy dominating the exchange 
coupling. At zero temperature, DTN undergoes a 
phase transition at a critical field Hi ~ 2.1 T above 
which nonzero spontaneous magnetization develops in 
the ground state and the spectrum of magnetic excita- 
tions becomes gapless. The system may then be thought 
of as a spin fluid formally described as a gas of hard- 
core bosons and the field-induced transition at Hi corre- 
sponds to Bose-Einstein condensation. A further transi- 
tion occurs at a second critical field H 2 ~ 12.6 T above 
which the ground state is a fully ordered ferromagnetic 
state. The field-induced quantum phase transitions de- 
scribed above have already attracted considerable ex- 
perimental interest through standard magnetometry, in- 
elastic neutron scattering, specific heat measurements, 
etc., followed by theoretical calculations based on suit- 
able Heisenberg models 1 . 

In particular, Zvyagin et alP^ have carried out de- 
tailed ESR measurements of magnetic excitations in a 
wide field range up to 25 T which includes the critical 
fields Hi and H 2 . The resulting picture was found to be 
generally consistent with early theoretical predictions^. 
In short, in the low-field region H < Hi the ground state 
carries zero azimuthal spin (S z = 0) while magnon exci- 
tations with S z = ±1 are separated by a gap which is un- 
ambiguously observed in the ESR spectrum through the 
expected AS Z = 1 transitions. In the high-field region, 
H > H2, the ESR spectrum is dominated at low tem- 
perature by AS Z = 1 transitions between a fully ordered 
ferromagnetic ground state and magnons that again ac- 
quire a nonzero energy gap. At finite temperature, a new 



feature appears in the ESR spectrum and corresponds to 
a AS Z — 1 transition between a magnon and a single-ion 
two-magnon bound state, a fact anticipated theoretically 
sometime agcBl A direct AS Z = 2 transition between the 
ordered ground state and a single-ion bound state is also 
observed. The physical picture becomes more involved 
in the intermediate field region Hi < H < H 2 but some 
progress has already been reported in recent literature^'. 

Our current task is to clarify some remaining dis- 
crepancies between theory and experiment. Thus 
we have carried out afresh a new set of ESR ex- 
periments performed in a frequency range 50-700 
GHz using a tunable-frequency submillimeter-wave ESR 
spectrometer^ equipped with Backward Wave Oscilla- 
tors as radiation sources and a 25 T resistive magnet. 
A transmission-type probe with a sample in the Fara- 
day geometry was employed (with the light propagation 
vector directed along the applied magnetic field H and 
the tetragonal c axis of the sample). High-quality sin- 
gle crystals of DTN with a typical size of about 3x3x5 
mm 3 (from a new batch, grown from aqueous solutions of 
thiourea and nickel chloride) were used. A silicon-based 
Dow Corning High Vacuum Grease 976 V was used to 
fix samples inside the probe. Particular attention was 
paid to measuring the temperature dependence of the ob- 
served ESR modes, especially in order to unambiguously 
resolve the contribution of two-magnon bound states in 
the high-field region H > H 2 . 

The main body of the paper is devoted to a theoretical 
analysis carried out in two steps. In Sec. II we adopt 
a quasi-one-dimensional Heisenberg model to calculate 
important features of the ESR spectrum through a sys- 
tematic strong-coupling expansion carried to third order. 
The general structure of the calculated spectrum agrees 
with experiment but important information concerning 
the relative intensity of the observed modes is practi- 
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FIG. 1: (color online). Frequency-field dependence of mag- 
netic excitations in DTN, with a uniform magnetic field H 
applied along the tetragonal c axis. Blue symbols denote ex- 
perimental data taken at T=1.7 K and red symbols at T=4.3 
K. Note that the mode E was observed in the Voight config- 
uration (with the light propagation vec tor directed perpen- 
dicular to the applied magnetic fleldjp^ while the rest of the 
modes were observed in the Voight as well as in the Faraday 
geometry. Solid lines correspond to results of calculations 
presented in Sec. II and are deliberately continued as dashed 
lines into the intermediate region Hi < H < H 2 - The loca- 
tion of critical fields Hi = 2.1 T and H2 = 12.6 T is indicated 
by vertical dashed lines. 



cally impossible to obtain within this essentially three- 
dimensional (3D) model. Thus, in Sec. Ill, we carry out 
such calculations within a strictly one-dimensional (ID) 
model through exact diagonalization on finite chains and 
a corresponding simulation of the relevant dynamic sus- 
ceptibilities. We are then able to analyze important fea- 
tures of the observed ESR spectrum over a wide field 
range including the intermediate region H\ < H < i?2- 
Our main conclusions arc summarized in Sec. IV. 



II. THREE-DIMENSIONAL MODEL 

The essential features of the observed ESR spectrum 
are illustrated in Fig. [l] together with some theoreti- 
cal predicti ons derived from a spin S = 1 Heisenberg 
HamiltoniarpJ^I 

H = J2 J »(Sr S 4+e J + £[£>(Sf ) 2 + g(i B HSn (1) 



where i denotes a generic site of a 3D lattice and e v 
with v = {x, y, z} count nearest neighbors. The exchange 
constants J v — {J x ,Jy,Jz\ may depend on the spe- 
cific lattice direction and are assumed to be significantly 
smaller than the easy-plane anisotropy (J u « D). Ac- 
tually, DTN is thought to be described by the quasi-one- 
dimensional limit of Eq. (TT~|) defined from J x = J y <C 



FIG. 2: (color online). A schematic view of the energy- 
momentum dispersions of magnetic excitations in an S — 1 
Heisenberg chain with strong easy-plane (D > 0) anisotropy 
for two typical fields H < Hi (left) and H > H 2 (right). Note 
that the ESR transitions denoted by A, B, C, E and F occur 
at k — 0, whereas transition G occurs at k = ir. Two-particle 
continua are not shown for simplicity. 



J z <C D, but the required theoretical analysis is essen- 
tially three-dimensional. Finally, an external magnetic 
field with strength H is applied in a direction perpendic- 
ular to the easy-plane. 

At zero field (H = 0) the ground state carries zero az- 
imuthal spin (S z — 0) and the magnon spectrum consists 
of two degenerate branches with S z = ±1 and energy- 
momentum dispersion uj = w(k) calculated through a 
systematic 1/D expansiorP carried to third order: 



w(k) = D + 2 J u cosk u 

V 

4 + J^cosK) 



1 

D 



1 



For nonzero but sufficiently low fields the S z = 
ground state remains unaffected while the degeneracy of 
the S z — ±1 magnon states is lifted (Fig. [2j left) to yield 
a twofold dispersion: 



w±(k) = w(k) ± g^BH 



(3) 



The ESR spectrum consists of two branches correspond- 
ing to A.S Z = ±1 transitions between the ground state 
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and k = magnons (modes A and B in Fig. [IJ. Thus 
the observed resonance frequencies are predicted to be: 

lo a = u> + gfJ, B H, lu b =uj - gfx B H (4) 

where ujq = uj(k = 0) is calculated from Eq. ([2|. Also 
note that the dispersion ui(k) of Eq. ^ exhibits a 
nonzero gap throughout the Brillouin zone, the smallest 
gap occurring at k = (tt,tt,tt). Therefore, the magnon 
frequencies of Eq. (|3| remain positive throughout the 
zone as long as H < Hi where Hi is a critical field de- 
fined from: 

gp B Hi = A , A = uj[k = (tt, tt, tt)] , (5) 

where the smallest gap A is again calculated from Eq. 
([2]) now applied for k = (71", 7T, tt). A corollary of the 
preceding discussion is that A < luq. 

When the field H exceeds its critical value Hi level 
crossing occurs and the azimuthal spin of the ground 
state no longer vanishes but increases with increasing 
field. Thus the system enters an intermediate phase 
through a field-induced quantum phase transition. The 
magnon spectrum is expected to be gapless in the inter- 
mediate phase but its detailed structure is now difficult 
to calculate. A systematic 1/D expansion is not feasible 
while semiclassical methods are generally inaccurate at 
strong anisotropy. Hence we postpone further discussion 
of the intermediate phase until Sec. Ill where a numerical 
calculation is carried out within a strictly ID model. 

The theoretical model of Eq. ([I]) becomes again 
tractable for sufficiently strong fields where the ground 
state is a completely ordered ferromagnetic state (Fig. 
|J right). The energy-momentum dispersion of single- 
magnon states is then given by 

e(k) = g» B H -D-2(J X + J y + J z ) 

+ 2(J x cosk x + JyCosky + J z cosk z ) (6) 

The lowest gap of this dispersion occurs at k = (tt, tt, tt) 
and is equal to gfi B H — D — A(J X + J y + J z ). Thus the 
ordered state is stable when the field exceeds a critical 
value given by 

g^ B H 2 = D + 4{J X + J y + J,) (7) 

Therefore, for H > H 2 , the ESR spectrum should be 
dominated by AS Z = 1 transitions between the com- 
pletely ordered ferromagnetic state and k = magnons. 
The resonance frequency is then calculated from 

lu c = e(k = 0) = gfi B H - D (8) 

and is found to be independent of the exchange constants. 

The physical picture is actually more involved for H > 
H 2 thanks to the appearance of an interesting class of 
two-magnon bound states. An exact calculation of such 
states is possible in the ID model through an elementary 
Bethe AnsatP. The two -magnon spectrum contains a 
"single-ion bound state" whose energy-momentum dis- 
persion extends well above the two-magnon continuum 



and was argued to be relevant for the analysis of the 
ESR spectrum observed in large- D systems^. 

However, a Bethe Ansatz is not applicable to the 3D 
model studied in this section. Thus we resort to a more 
direct method developed long time ago by Wortis^l for 
the calculation of two-magnon bound states in ferromag- 
nets with arbitrary lattice dimension. The method is here 
generalized to account for easy-plane anisotropy with 
strength D and is employed in conjunction with the 1/D 
expansion when analytical treatment is no longer feasible. 
Thus we were able to calculate the energy-momentum 
dispersion of the single-ion mode E — E(k) to third or- 
der in the 1/D expansion. 

We defer for the moment discussion of a k = 0, 
AS Z — 2 transition between the ordered ground state 
and a single-ion bound state. Instead, we turn our atten- 
tion to AS Z = 1 transitions between single magnons and 
single-ion bound states. These are absent at zero temper- 
ature but may occur with nonvanishing intensity at finite 
temperature. The corresponding resonance frequencies 
are then given by w(k) = E(k) — e(k) where k extends 
over the entire Brillouin zone. Hence, at sufficiently low 
but nonzero temperature, resonance frequencies are ex- 
pected to be observed throughout a band top < lo < ljq 
where the lower frequency is calculated to third order: 

lu f = E(k = 0) - e(k = 0) 



= gi x B H + D-A{J x + J y + J z ) 




whereas the upper frequency is given by 

loq = E[k = (n, 7r, 7r)] — e[k = (tt, tt, tt)] 

= g f i B H + D (10) 

which is an exact result independent of the exchange con- 
stants, in analogy with the resonance frequency ljc of Eq. 

To summarize, the single-ion (FG) band is absent at 
zero temperature but acquires nonvanishing intensity at 
finite temperature. As was argued in Ref. 0] and further 
discussed in Sec. Ill of the present paper, the inten- 
sity is expected to display a characteristic double peak 
as a function of frequency at fixed external field, or as 
a function of field at fixed frequency. Therefore, both 
frequencies uif and loq are associated with the single- ion 
bound state and are relevant for the analysis of actual 
experiments. In this respect, it is worth mentioning that 
the absorption corresponding to mode G was initially ob- 
served in previous experiments^. This absorption was in- 
terpreted as an artifact originating in the superficial layer 
of DTN crystals attacked by a GE- varnish solvent used 
to fix the sample within the sample holder. 

However, our theoretical analysis suggests that the G 
mode is actually an inseparable partner in a doubly- 
peaked FG band associated with the single-ion bound 
state. Indeed, our current experiment supports such an 
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FIG. 3: (color online). ESR transmittance spectra in DTN 
taken at frequency 647 GHz for three representative temper- 
ature values. Vertical line C indicates the location of the 
calculated single-magnon resonance, while F and G are the 
calculated boundaries of the single-ion (FG) band. Note that 
experiments were performed in magnetic fields up to 25 T 
and thus the single-magnon (C) resonance is not shown in 
this figure. 

interpretation, as shown in Fig. [3] where the transmit- 
tance measured at fixed frequency / = 647 GHz displays 
a characteristic double peak as a function of the applied 
field. Also note that the double peak is uneven with most 
power absorbed for frequencies near the G boundary, an 
experimental fact that will be shown to be consistent 
with a numerical calculation of power absorption in Sec. 
III. Another important feature of Fig. [3] is the appar- 
ent vanishing of intensity at relatively low temperatures 
(e.g., T = 1.5 K), an experimental fact that is consis- 
tent with our interpretation of the FG resonance band as 
the result of transitions between excited states; namely, 
transitions between single magnons and single-ion two- 
magnon bound states (Fig. [2j right). The overall picture 
suggested by Fig. [3] is robust in the high-field region 
(H > H 2 ) and extends into the intermediate and low- 
field region (H < H 2 ) as shown in Fig. [4] and further 
discussed in Sec. III. 

Finally, we return briefly to the possibility of a AS Z = 
2 transition between the ordered ground state and a k = 
single-ion bound state, which would lead to a resonance 
frequency 

uj e = E(k = 0) = uj f + uj c (11) 

but zero intensity thanks to the axial symmetry adopted 
in our theoretical models. However, crystal symmetry is 
compatible with some deviations from strict axial sym- 
metry which apart from a tiny field misalignment may 
render mode E observable. In fact, such a mode was 
previous ly o bserved in DTN with a sample in the Voigt 
geometry^ and is included in Fig. [lj 

The remainder of this section is devoted to a brief 
discussion concerning the choice of suitable parameters. 
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FIG. 4: (color online). ESR transmittance spectra for three 
characteristic frequencies and a wide field range up to 25 T. 
Note that the specific choice of frequencies is such that all 
possible modes appear in the figure. 

The simplest possibility is to fit the zero-field magnon 
dispersion given by Eq. ^ to the dispersion measured 
via inelastic neutron scattering^. A good fit is obtained, 
see Fig. [5| with the choice of parameters^ 

D = 7.72 K, J x = J y = 0.2 K, J z = 1.86 K (12) 

In particular, ujq = uj(k = 0) = 257 GHz, in fair agree- 
ment with the experimental value ujq = 267 GHz inde- 
pendently obtained through ESR. Consequently, the the- 
oretically predicted branches A and B in the ESR spec- 
trum agree with experiment if we further choose a gyro- 
magnetic ration g = 2.22. The critical field calculated 
from Eq. (J5j) is Hi = 2.08 T, in excellent agreement with 
the experimental value H± =2.1 T. 

Thus the preceding choice of parameters yields a suf- 
ficiently accurate description of the low-field region H < 
Hi . But such a choice leads to poor quantitative predic- 
tions in the high- field region H > H 2 . For example, the 
critical field H 2 calculated from Eq. Q is H 2 = 11.24 
T, to be compared with the experimental H 2 — 12.6 T. 
Similarly, the exact magnon branch u>c of Eq. (fsl) sub- 
stantially disagrees with experiment when D = 7.72 K. 

Instead, an excellent fit of mode C is obtained using 
the parameters^ 

3 = 2.22, L> = 8.9K (13) 

We adopt these values and fix the remaining (exchange) 
constants via a least-square fit of the zero-field dispersion 
of Eq. ([2| to the experimental dispersion, see Fig. [5| to 
obtain 

J x = J y = 0.34 K, J z = 1.82 K (14) 

which arc significantly different from J x = J y = 0.18 K 
and J z = 2.2 K obtained in Ref. fusing a self-consistent 
semiclassical method to calculate the zero-field magnon 
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FIG. 5: DTN dispersion of magnetic excitations at zero field 
calculated along three directions in the Brillouin zone using 
Eq. Q a nd parameters taken from Eq. ( 12 1 (dashed lines) 
and Eqs. ||13[)-(14| (solid lines). Symbols denote inelastic neu- 



tron scattering data taken from Ref. [T] Energy is measured 
in degrees K. 



dispersion. Here, to be consistent, we employ the pa- 
rameters of Eq. (131 and Eq. (14 1 to calculate the crit- 
ical fields Hi = 2.30 T and H 2 = 12.68 T which are in 
rough agreement with the experimental Hi = 2.1 T and 
Hi = 12.6 T. We also employ the same parameters to 



calculate the branches of the ESR spectrum shown by 
straight lines A, B, C, F, G, E in Fig. [IJ which are again 
in rough agreement with experiment. A schematic repre- 
sentation of the ESR transitions that correspond to the 
above modes is given in Fig. [2] 

In our opinion, a completely satisfactory quantitative 
agreement with experiment should not be expected for 
two reasons. First, crystal symmetry of DTN allows, 
in principle, inclusion of several more terms of unknown 
strength in the effective Heisenberg Hamiltonian of Eq. 
such as Dzyaloshinskii-Moriya interactions between 
neighboring Ni sites 5 . Second, even within the limits of 
model of Eq. , accurate quantitative predictions may 
be difficult to obtain within the limited third-order 1 JD 
expansion, while similar uncertainties may arise in the 
context of sophisticated semiclassical methods^. Never- 
theless, we believe that the overall qualitative picture is 
substantially correct. 



III. ONE-DIMENSIONAL MODEL 

Yet there are several questions that are difficult to set- 
tle within the essentially 3D model of Eq. ([I]), such as 
the calculation of intensities of the various ESR modes, 
the structure of the spectrum in the intermediate phase, 
etc. But the essential features of the observed spectrum 
are already accounted for by a relatively simple ID model 
which is further employed in the present paper to discuss 
the remaining questions. 

In order to establish consistency with the earlier worlP 
we first restrict the main results of Sec. II to the ID 
model through the formal substitution 



J x Jy 



0, 



(15) 



For instance, the Heisenberg Hamiltonian of Eq. ([IJ now 
reads 



N 



H = ]T[J(S„ • S n+1 ) + D(S z n f + gnaHSZ] (16) 

71=1 

where the summation extends over the N sites of a ID 
lattice assumed to be periodic. 

The zero-field magnon dispersion of Eq. ([2| reduces to 



w(Jfe) 



D 2 



(1 + 2sin 2 k) 



D 1 + —cosh 

J3 i 

-jy^[2sin 2 k — — (1 + 8sin 2 k)cosk] 



(17) 



and agrees with an early calculation within the ID 
modePM Recently, several mor e term s beyond the 
third order have become available 1 12 * 13 !. However, for 
anisotropy strengths of current interest, the third-order 
result (17) proves to be sufficiently accurate. 



In particular, the ESR modes A and B may be calcu- 
lated from Eq. Q now applied with 



wo = u{k = 0) = D 1 + — + 



2.J 
~D 



J2_ 
D 2 



J 3 
2D 3 



(18) 



Similarly, the critical field Hi is given by Eq. |5]) applied 
with 



A = u>(k = 7r) = D 1 - 



2J 

15 



D 2 



J 3 
2D 3 



(19) 



Analogous results may be obtained for sufficiently strong 
fields in the region H > H 2 where the ground state is 
a fully ordered ferromagnetic state. The single-magnon 
dispersion of Eq. ^ reduces to 



e(k) = g\x B H - D- 2J(1 - cosk) 



(20) 



whose lowest gap occurs at k — tt and is equal to gfi B H - 
D — 4 J. Therefore, the upper critical field is given by 



g^ B H 2 =D + 4J, 



(21) 



which agrees with the ID reduction of Eq. ([T]). For 
H > H 2 , the domimant resonance frequency ljc arises 
from AS Z — 1 transitions between the ordered state and 
k = magnons: 

uj c = e{k = 0)=gfi B H-D (22) 

which is exact and coincides with the 3D result of Eq. 

The dispersion of the single-ion two-magnon bo und 
state can be calculated exactly within the ID modeP^ 
but a third-order approximation is sufficient for our pur- 
poses: 



E(k) = 2gn B H + D 



_4J 2J2 



D 



J > n k 

D~ 3)C0S 2 



(23) 
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Thus a AS Z = 1 transition between a k — magnon 
and a k = single-ion two-magnon bound state yields a 
resonance frequency 



Intensity I(f,h,x) 



uj f = E(k = 0) - e(fc 



1 - 



0) 
4J 
~D 



2J2 



D 3 



(24) 



which agrees with the ID reduction of Eq. ([9| , whereas 
a transition between a k = n magnon and a k = 7r single- 
ion bound state yields 



ujq = E(k = 7r) — e(fc = 7r) = q^bH + D 



(25) 



which is exact (independent of J) and coincides with the 
3D result of Eq. ( 10 ). Similar transitions occur for other 



values of k throughout the Brillouin zone and lead to 
a band of resonance frequencies in the region ujp < uj < 
ujq- Although the intensity of such transitions vanishes at 
zero temperature, nonzero intensity is expected to occur 
at finite temperature, an issue to be discussed in detail 
in the continuation of this section. 

First, a digression concerning the choice of parame- 
ters within the ID model. Recall that the single-magnon 
resonance frequency ujc given in Eq. ( 22 ) is an exact pre- 
diction of the ID as well as the 3D model; see Eq. ([8]). 
Therefore, we adopt in this section the choice of the gy- 
romagnetic ratio g and anisotropy D already made in Eq. 
( 13 1 and the only remaining parameter is the exchange 



constant J or, equivalently, the dimensionless ratio J/D. 
A semi-quantitative agreement with experiment is ob- 
tained with the choice 



2.22, 



D = 8.9 K, 



J 
D 



(26) 



which will be adopted in all calculations presented in this 
section. For convenience, we use rationalized variables 
such that frequency / — lo/2-k is measured in units of 
D/2-kH = 185.45 GHz, magnetic field h = g^ B H/D in 
units of D/gfiB = 6 T and temperature r = T/D in 
units of D = 8.9 K. 

The main issue addressed in this section is an explicit 
calculation of power absorption. In a typical ESR exper- 
iment a microwave field of angular frequency co is applied 
in the basal plane along, say, the x-axis, in addition to 
a uniform bias field H applied along the z-axis. The in- 
tensity or power absorption per site is defined up to an 
overall multiplicative constant by 



(27) 



where the imaginary part of the susceptibility is given 
by23 



a , b 



(28) 

Here sums extend over all eigenstates \a) of Hamiltonian 
of Eq. (16), E a are the corresponding eigenvalues, j3 = 



Magnetic Field h 

FIG. 6: (color online). Colored surface represents the normal- 
ized ESR intensity /(/, h, r) calculated for a spin chain with 
N=10 as a function of frequency / and magnetic field h, at 
fixed temperature r = 0.2 (T — 1.8 K). Solid lines corre- 
spond to results of calculation within the ID model of Sec. 
Ill and are deliberately extended as dashed lines into the in- 
termediate region hi < h < h%. The location of critical fields 
hi — 0.57 and h% = 2 is indicated by vertical dashed lines. 
Note that /, h and r are all measured in rationalized units 
defined after Eq. pHl. 



1/T is the inverse temperature, and Z = ^2 a e~^ Ea is 
the total partition function. Finally, matrix elements in 
Eq. ( |28[ ) involve the total spin operator in the x-direction 
/; — "V ^ 

An analytical calculation of x"( w ) is on f °f question 
except in very special limits^. We thus resort to a nu- 
merical calculation based on Eq. (28) and a complete 



diagonalization of the spin-1 Hamiltonian of Eq. (16) 



defined on a finite periodic chains with size N as large 
as 12. Actually, explicit results presented below were 
obtained on a chain with N=10, whereas N=12 chains 
were occasionally used for consistency checks. On a fi- 
nite chain Eq. ( 28 ) yields a susceptibility that is a sum 
of weighted (5-functions and is thus rather spiky. Hence 
we adopted an empirical smoothing process to obtain an 
intensity 



I = I(f,h,r) 



(29) 



that is a reasonably smooth function of frequency /, mag- 
netic field h, and temperature r, measured in rationalized 
physical units defined in the text following Eq. ( 26 1 . Our 



main task is then to analyze the calculated intensity as 
a function of all three variables. 

In Fig. m we present our results for the (colored) sur- 
face I = l(f, h) at fixed temperature r — 0.2 (T = 1.8 
K) which is a typical relatively low temperature of ex- 
perimental interestP. Superimposed in the same figure 
are the analytical predictions for the two critical fields 
hi = 0.57 and hi — 2 calculated from Eq. (19) and Eq. 



(21 ) adapted to rationalized units, as well as correspond- 
ing predictions for the resonance lines A, B for h < hi 
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Intensity I(f,h,t) 



Intensity I(f,h,T) 



Magnetic Field h 

FIG. 7: (color online). Same as Fig. [6]but now focusing on 
the low-frequency end of the intermediate phase h\ < h < hi . 
Note the formation of a V-like structure with rapidly decreas- 
ing intensity as one approaches the center of the intermediate 
phase. 



and C, F, G for h > hi calculated earlier in this section. 
Note that mode E (see Fig. [T]) does not appear in Fig. 
[6] because it corresponds to a AS Z = 2 transition and its 
intensity / vanishes within the strictly axially symmetric 
model of Eq. (16). Several important facts have already 



become apparent in Fig. [6] which we analyse in turn: 

(a) We note that the magnon resonance lines A, B and C 
roughly coincide with the maxima of the calculated 
intensity as expected in the low-temperature region. 
Nevertheless, the chosen temperature T = 1.8 K is 
sufficiently high to account for the anticipated line 
broadening which is also apparent in Fig. [6j Yet, this 
temperature is too low to yield a significant signal for 
the single-ion bound state, as shown in Fig. [6] where 
the intensity practically vanishes in the FG region. 

(b) We examine the results of Fig. [6] in the intermediate 
region h\ < h < hi where analytical predictions are 
practically absent. The most conspicuous feature is a 
tail of line G with strong intensity in the intermediate 
region, which persists even at very low temperature 
where line G itself looses its intensity for h > hi. 
Therefore, the G-tail corresponds to some sort of a 
collective excitation that is robustly present in our 
current experiment and requires further theoretical 
investigation. On the other hand, mode F acquires a 
tail into the intermediate region with intensity that 
diminishes at low temperature and is thus invisible 
in Fig. [6j 

(c) We note that the magnon lines A, B and C also 
acquire tails but with intensity that gradually van- 
ishes as one approaches the center of the intermediate 
phase. The structure of the tails becomes apparent in 
Fig. [7] which focuses on the low- frequency end of the 
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FIG. 8: (color online). Same as Fig. |6]but with intensity now 
calculated at a much higher temperature r = 2 (T = 18 K). 
Note that we concentrate on the high-field region h > hi = 2 
(H > 12 T) in order to emphasize the significant enhancement 
of intensity in the FG band (especially near the G boundary) 
which provides unambiguous evidence for the existence of a 
single-ion two-magnon bound state. 



intermediate phase. Thus we reveal a V-like struc- 
ture with intensity that gradually vanishes as one ap- 
proaches the center. This picture apparently contra- 
dicts the result of Cox et alP who predict by a similar 
calculation a Y-like structure with intensity that re- 
mains finite and practically constant near the center. 
On the other hand, our result is consistent with a 
rounding of a V into a U structure predicted to oc- 
cur in the presence of a small Dzyaloshinskii-Moriya 
anisotropy treated by a semiclassical methocP. 

(d) As mentioned already, Fig. [6] as well as experi- 
ment indicate absence of measurable intensity in the 
single-ion (FG) band at the relatively low tempera- 
ture r = 0.2 (T = 1.8 K). However, the FG band 
is significantly activated at higher temperature, as 
demonstrated in Fig. [8] which displays the intensity 
in the field region h > hi = 2 (H > 12 K) and tem- 
perature t = 2 (T = 18 K). The calculated FG band 
is highly populated at this temperature with most of 
the intensity concentrated near the G boundary. 

To understand the preceding result in some detail, we 
depict in Fig. [9] the intensity as a function of frequency 
at a fixed field h = 2.5 (H = 15 T) and selected val- 
ues of temperature. At the lowest temperature r = 0.1 
(T = 0.9 K) considered in Fig. [9j the dominant feature 
is the magnon resonance C while there is no sign for a 
single-ion bound state. On the contrary, the FG band 
is activated already at temperature r = 0.5 (T = 4.5 
K) employed in actual experiments^^. The FG signal is 
further enhanced at higher temperature, as is evident in 
the t = 1 (T = 9 K) and t — 2 (T — 18 K) entries. 
Also evident is the formation of a double peak in the 
FG region, with the dominant peak occurring near the G 
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FIG. 9: Calculated normalized intensity /(/, h, r) scaled with 
temperature r, as a function of frequency / at a fixed field 
h — 2.5 and four typical values of temperature r. Vertical 
lines C, F and G indicate the location of the single-magnon 
resonance C and the boundaries of the single-ion (FG) band 
calculated within the ID model. Recall that f,h and r are all 
measured in rationalized units defined after Eq. ( 26 1. 



FIG. 10: Total (integrated) intensity of the FG band as a 
function of temperature at fixed frequency / = 647 GHz cal- 
culated within the ID model. Symbols denote experimental 
data extracted from field integration of ESR spectra. 



interpretation of the ESR signal of the single-ion two- 
magnon bound state. 



boundary while a peak of lower intensity develops near 
the F boundary. The relative enhancement of the inten- 
sity near the G boundary is likely due to the fact that it 
involves transitions between k = tt single magnons and 
k = tt single-ion bound states, where the magnon ac- 
quires its lowest gap and is thus more heavily populated 
at finite temperature than, say, k — magnons. 

The preceding theoretical findings are consistent with 
experimental results of the type shown in Fig. [3j with 
due attention to the fact that Fig. [3] depicts the trans- 
mittance as a function of applied field at fixed frequency. 
In any case, both theory and experiment suggest a dom- 
inant peak near the G boundary followed by a secondary 
peak (a knee) near the F boundary. The two peaks are 
partners in a doubly-peaked FG band that cannot be 
separated in any meaningful way. Thus it is difficult to 
measure or calculate their relative intensity. Neverthe- 
less, it is possible to calculate the total intensity of the 
FG band: 



'FG 



J a 



h, r)dh 



(30) 



where integration extends over a field interval a < h < b 
chosen empirically so that it encompass the entire FG 
band. The total intensity of Eq. (30) is depicted in 
Fig. [lO] as a function of temperature at fixed frequency 
/ = 647 GHz, together with experimental results ob- 
tained by applying a similar integration process to data 
of the type shown in Fig. [3j Taking into account that 
intensity is displayed in "arbitrary units" , the qualitative 
agreement between theory and experiment shown in Fig. 
10 is satisfactory and fully consistent with our current 



IV. CONCLUSIONS 

As far as the general structure of the observed ESR 
spectrum is concerned, the theoretical predictions of the 
3D model of Eq. Q and the ID model of Eq. Q 
are qualitatively similar and in reasonable quantitative 
agreement with experiment. But a detailed investigation 
of the remaining discrepancies required a calculation of 
the intensities of the various ESR modes, which is not 
feasible within the 3D model. Thus most of our effort 
was devoted to a detailed numerical calculation of inten- 
sity within the ID model. The main new results are the 
following: 

While there have been numerous theoretical predic- 
tions for the occurrence of two-magnon bound states 
in quantum spin systems, experimental observation has 
been rather slow. Perhaps, the most interesting feature 
of the ESR spectrum in large- D systems is the evidence 
it provides for the existence of the so-called single-ion 
two-magnon bound states. The original theoretical sug- 
gestion was made some time agc^ and was thought to 
explain ESR data obtained on a large- D compound ab- 
breviated as NENO^K^. But a thorough experimental 
investigation was carried o ut m ore recently in relation to 
the title compound (DTN^pH. 

Our present investigation clearly suggests that the F 
and G lines are inseparable partners in a doubly-peaked 
FG band which originates in transitions between single 
magnons and single-ion two-magnon bound states. In 
fact, the G mode absorbs most of the intensity and is thus 
far from extraneous. This mode is especially interesting 
in relation to the fact that the corresponding resonance 
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line ojq = g/isH + D is an exact theoretical prediction 
both within the 3D model of Eq. and the ID model 
of Eq. fl6l); see Eqs. (pi) and pHl. 



Our numerical calculation also sheds light on the struc- 
ture of the magnetic excitation spectrum in the interme- 
diate phase Hi < H < H2 where analytical results are 
practically absent. As is evident in Fig. [6j a tail of line G 
with strong intensity survives in the intermediate region 
even at low temperature where line G itself looses its in- 
tensity for H > H2- Such a tail should thus be attributed 
to a high-frequency collective excitation that appears in 
the intermediate phase as a shadow of the single-ion two- 
magnon bound state, an issue that deserves further the- 
oretical attention. 

Finally, the current calculation does not support the 
occurrence of a low-frequency Y structure suggested by 
Cox et alP, even though they also employ the ID model 
of Eq. (16 1 to calculate the susceptibility x"(w). In fact, 
we find a V structure with rapidly decreasing intensity 
near the center of the intermediate phase. As such the V 
structure is expected to be especially vulnerable to small 



perturbations that are ever present in effective Heisen- 
berg models. This may explain the deformation of the V 
into a U shape in the presence of a small Dzyaloshinskii- 
Moriya anisotropj0. 



ACKNOWLEDGEMENTS 

The theoretical part of this work has been supported 
by the INT 238475 LOTHERM project. The experimen- 
tal part was performed at the National High Magnetic 
Field Laboratory, Tallahassee, FL, which is supported 
by NSF Cooperative Agreement No. DMR-0654118, by 
the State of Florida, and by the DOE. S.Z. appreciates 
the support of the Deutsche Forschungsgemeinschaft and 
EuroMagNET II (EU contract No. 228043). We are 
grateful to Dr A.V Sizanov and Dr J. Oitmaa for useful 
correspondence concerning the results of Refs. l51andH"2l 
respectively. 



1 V. S. Zapf, D. Zocco, B. R. Hansen, M. Jaime, N. Harrison, 
C. D. Batista, M. Kenzelmann, C. Niedermayer, A. Lac- 
erda, and A. Paduan-Filho, Phys. Rev. Lett 96, 077204 
(2006) 

2 S. A. Zvyagin, J. Wosnitza, CD. Batista, M. Tsukamoto, 
N. Kawashima, J. Krzystek, V. S. Zapf, M. Jaime, N. 

F. Oliveira, Jr., and A. Paduan-Filho, Phys. Rev. Lett 98, 
047205 (2007) 

3 S. A. Zvyagin, C. D. Batista, J. Krzystek, V. S. Zapf, 
M. Jaime, A. Paduan-Filho, and J. Wosnitza , Physica 
B 403, 1497 (2008) 

4 N. Papanicolaou, A. Orendacova, and M. Orendac, Phys. 
Rev. B 56, 8786 (1997) 

5 S. A. Zvyagin, J. Wosnitza, A. K. Kolezhuk, V. S. Zapf, 
M. Jaime, A. Paduan-Filho, V. N. Glazkov, S. S. Sosin, 
and A. I. Smirnov, Phys. Rev. B 77, 092413 (2008) 

6 S. Cox, R. D. McDonald, M. Armanious, P. Sengupta, and 
A. Paduan-Filho, Phys. Rev. Lett. 101, 087602 (2008) 

7 S. A. Zvyagin, J. Krzystek, P.H.M. van Loosderecht, 

G. Dhalenne, and A. Revcolewschi, Physica B 346-347, 
1 (2004) 



8 A. V. Sizanov and A. V. Syromyatnikov, 
arXiv: 1 105 . 1077vl [cond-mat .str-el] 

9 N. Papanicolaou and G. C. Psaltakis, Phys. Rev. B 35, 
342 (1987) 

10 M. Wortis, Phys. Rev. 132, 85 (1963) 

11 N. Papanicolaou and P. N. Spathis, Phys. Rev. B 52, 16001 
(1995) 

12 A. F. Albuquerque, C. J. Hamer, and J. Oitmaa, Phys. 
Rev. B 79, 054412 (2009) 

13 J. Oitmaa, private communication. 

14 C. P. Slichter, Principles of Magnetic Resonance (Springer- 
Verlag, Berlin, 1978) 

15 S. A. Zvyagin, V. V. Eremenco, V. V. Pishko, A. Feher, 
M. Orendac, and A. Orendacova, Low Temp. Phys. 21, 
680 (1995) 

16 S. A. Zvyagin, T. Rieth, M. Sieling, S. Schmidt, and 
B. Luthi, Czech. J. Phys. 46, 1937 (1996) 

17 A. K. Kolezhuk, and H.-J. Mikeska, Phys. Rev. B 65, 
014413 (2001) 



